DESCRIPTION OF THE SCENARIO MACHINE 
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ABSTRACT 

We present here an updated description of the "Scenario Machine" code. This tool is 
used to carry out a population synthesis of binary stars. Previous version of the descrip- 
tion can be found at |http://xray.sai.msu.ru^mystery//articles/review/contents.html[ see also 



(jLipunov et al. Hl996lM 7 



Subject 



binaries: close — binaries: general 



1. Basic equations and initial distributions 

We use the current scenario of evolution of bi- 
nary stellar systems based up on the original idea s 
that appeared in the papers b y lPaczvn'ski ( 1971 ): 



Tutukov fe Yungelson I d 19731) 



van den Heuvel fc Heise ( 19721) (see also re- 
view by Ivan den Heuvel IS)). The sce- 



nario for normal star evolution was joined with 
the ideas of ne utron star e v olutio n (see pi- 
oneer works by Shvartsman (Il970 , 1971allb ): 



Illarionov & Sunyaev (119751): IShakura 
Bisnovatvi-Kogan fc Romberg I (|l976l ), 
Lipunov fc Shakura I (1 19761) 



(|l975l ): 



Savoniie fc van den Heuvel (|l977l )V This joint 



scenario has allowed to construct a two-dimensional 
classification of possible states of bina r y syste ms 
containing NS ( Kornilov fc Lipunov ( 1983al lbh; 



Lipunov I (|l992l) ). According to this classification, 
we will distinguish four basic evolutionary stages 
for a normal star in a binary system: 



I — A main sequence (MS) star inside its Roche 
lobe (RL); 

II — A post-MS star inside its RL; 

III — A MS or post-MS star filling its RL; the mass 
is transferred onto the companion. 



1 E-mail 

2 E-mail 

3 E-mail 
4 E-mail 



lipunov@xray.sai.msu.ru 
pk@. sai.msu.ru 
mystery@xray.sai.msu.ru 
a78b@yandex.ru 



IV — A helium star left behind the mass-transfer 
in case II and III of binary evolution; may 
be in the form of a hot white dwarf (for M < 
2.5M Q ), or a non-degenerate helium star (a 
Wolf-Rayet-star in case of initial MS mass 
> 1OM ). 



The evolution of single stars can be represented 
as a chain of consecutive stages: I — >• II — >• compact 
remnant; the evolution of the most massive single 
stars probably looks like I — > II — » IV — * compact 
remnant. The component of a binary system can 
evolve like I — > II — > III — > IV — > compact remnant. 

In our calculations we choose the distributions 
of initial binary parameters: mass of the primary 
zero age main sequence component (ZAMS), Mi, 
the binary mass ratio, q = M2/M1 < 1, the orbital 
separation a. Zero initial eccentricity is assumed. 

The distribution of binaries by o rbital separa- 
tions can be take n from observations (jKraicheva et al. 
198lUAbtl[l98l . 



/(log a) = const, 

max(10i? Q ,RL[Afi]) < a < 10 6 i? G 



(1) 



Especially important from the evolutionary point 
of view is how differe nt are in i tial m asses of the 
components (see e.g. I Trimble I (|1983l) ). We have 
parametrized it by a power-law shape, assuming 
the primary mass to obey Salpeter's power law: 



f(M) 



M -2.36 j 



0.1M Q < Mi < 12OM , (2) 



1 



/(<?) 



q = Mi/Mi < 1; 



(3) 



We should note that some apparently reason- 
able distributions - such as both the primary and 
secondary mass obeying Salpeter's law, or "hier- 
archical" distributions involving the assumption 
that the total binary mass and primary's mass are 
distributed according to the Salpeter mass func- 
tion - all yield essentially flat-like distributions by 
the mass ratio (i.e. with our parameter a q ~ 0). 

We assume that the neutron star is formed 
in the core collapse of the pre-supernova star. 
Masses of the young neutron stars are randomly 
distributed in the range M™ s n - M%§*. Initial 
NS masses are taken to be in the range M^s — 
1.25 — 1.4M©. The range of initial masses of the 
young NSs was based on the masses of the neu- 
tron star in the B1913+16 binary system and of 
the radio pulsar in the J0737-3039 binary. In the 
B1913+16 system (the Hulse- Taylor pulsar, radio 
pulsar + neutron star) the mass of the neutron 
star, which is definitely not accreting matter from 
the optical donor, is M NR = 1 3873 ± 0.000671^ 



0: 



dThorsett fc Chakrabartv1ll999t lW ex et al. l20Qi 
Weisberg fc Taylor 1 12003|) . The mass of the pul- 
sar in the J0737-3039 system (radio pulsar + radio 
pulsar), which likewise does not accrete from an 
optical companion , is Mpsr = 1.250 ± O.OIOMq 



( Lyne et al. 2004 ). The mass range of stars pro- 



ducing neutron stars in the end of their evolution is 
assumed to be M n - M&; stars with initial masses 
M > Mb are assumed to leave behind black holes; 
M n is taken to be equal to 10M© in most cases 
(but in general this parameter is free). Note that 
according to some stellar models, a very massive 
star (« 50 — 100M Q ) can leave behind a neutron 
star as a remnant due to very strong mass loss 
via powerful stellar wind, so we account for this 
possibility in the corresponding models. 

We take into account that the collapse of mas- 
sive star into a neutron star can be asymmetrical, 
so that the newborn neutron star can acquire an 
additional, presumably randomly oriented in space 
kick velocity w (see Section |4] below for more de- 
tails). 

The magnetic field of rotating compact ob- 
jects (neutron stars and white dwarfs) largely 
define the evolutionary st age of the compac t 
object in a binary system ( Shvartsman I ( 1970l ): 



Davidson fe Ostriker I (|l973f ); IHlarionov fe Sunvaev 
(|l975l )). so we use the general classificati on of mag- 



netic rotating compact objects (see e.g. iLipunov 



(1992)) m our calculations. The initial magnetic 
dipole moment of the newborn neutron star is 
taken according to the distribution 



/(log/i) cx const, 

10 28 < n < 10 32 G cm 3 , 



(4) 



The initial rotational period of the newborn neu- 
tron star is assumed to be ~ 10 ms. 

It is not definitely clear as yet whether the mag- 
netic field of neutron s tars decays or not (se e for a 



comprehensive review IChanmugam I (|l992h 'l. We 



assume that the magnetic fields of neutron stars 
decays exponentially on a timescale of td (usually 
we take this parameter to be equal to 10 8 , 5 • 10 7 
and 10 7 years). A radio pulsar is assumed to be 
"switched on" until its period P (in seconds) has 
reached the "death-line" defined by the relation 
A*3o/^d = 0.4, where /X30 is the dipole magnetic 
moment in units of 10 30 G cm 3 . 

We assumed that magnetic fields of neu- 
tron stars decay exponentially to minimal value 
B m in — 8 • 10 7 G and do not decay further: 



B 



B exp(-i/i rf ), t < t d ln(B /B min ), 
B m i n ,t > td \n(B n /B min ). 



(5) 



Parameters Bq and td in equation ([5]) are the ini- 
tial field strength and the field decay time. 

We also assume that the mass limit for neutron 
stars (the Oppenheimer-Volkoff limit) is Mov — 
2.OM0 (in general, it is a free parameter in the 
code; it depends on equation of state of the mate- 
rial of the neutron star). 

The most massive stars are assumed to leave be- 
hind black holes after the collapse, provided that 
the progenitor mass before the collapse has a mass 
M cr . The masses of the black holes are calcu- 
lated as Mbh — kbhMp re sN, where the parameter 
kbh = 0.0 — 1.0, Mp re sN is the mass of the pre- 
supernova star. 

We consider binaries with Mi > O.8M with a 
constant chemical (solar) composition. The pro- 
cess of mass transfer between the components is 
treated as conservative when appropriate, that is 
the total angular momentum of the binary system 



2 



is assumed to be constant. If the accretion rate 
from one component to another is sufficiently high 
(say, the mass transfer occurs on a timescale few 
times shorter than the thermal Kelvin-Hclmholtz 
time for the normal companion) or a compact ob- 
ject is engulfed by a giant companion, the com- 
mon envelope stage of binary evolution can begin 
( Paczvnski 1 19761 : Ivan den Heuvel I 



1983). 



Other cases of non-conservative evolution (for 
example, stages with strong stellar wind or those 
where the loss of binary angular momentum oc- 
curs due to gravitational radiation or magnetic 
stellar wind) are trea ted using the well know n 
rescriptions (see e. g. IVerbunt fc Zwaanl (119 81): 



(1988)) 



Rap paport et al. I (|1982l ): iLipunov fc Postnov 



2. Evolutionary scenario for binary stars 

Significant discoveries in the X-ray astronomy 
made during the last decades stimulated the as- 
tronomers to search for particular evolutionary 
ways of obtaining each type of observational ap- 
pearance of white dwarfs, neutron stars and black 
holes, the vast majority of which harbours in bi- 
naries. Taken as a whole, these ways costitute a 
general evolutionary scheme, or the "evolutionary 
scenario". We follow the basic ideas about stel- 
lar evolution to describe evolution of binaries both 
with normal and compact companions. 

To avoid extensive numerical calculations in 
the statistical simulations, we treat the continu- 
ous evolution of each binary component as a se- 
quence of a finite number of basic evolutionary 
states (for example, main sequence, red super- 
giant, Wolf-Rayet star, hot white dwarf, etc.), at 
which stellar parameters significantly differ from 
each other. The evolutionary state of the binary 
can thus be determined as a combination of the 
states of each component, and is changed once the 
more rapidly evolving component goes to the next 
evolutionary stage. 

At each such stage, we assume that the star 
does not change its physical parameters (mass, ra- 
dius, luminosity, the rate and velocity of stellar 
wind, etc.) that have effect on the evolution of 
the companion (especially in the case of compact 
magnetized stars) . Every time the faster evolving 
component passes into the next stage, we recalcu- 
late its parameters. Depending on the evolution- 



ary stage, the state of the slower evolving star is 
changed or can remain unchanged. With some ex- 
ceptions (such as the common envelope stage and 
supernova explosion), states of both components 
cannot change simultaneously. Whenever possible 
we use analytical approximations for stellar pa- 
rameters. 

Prior to describing the basic evolutionary states 
of the normal component, we note that unlike 
single stars, the evolution of a binary compo- 
nent is not fully determined by the initial mass 
and chemical composition only. The primary star 
can fill its Roche lobe either when it is on the 
main sequence, or when it has a (degenerate) he- 
lium or carbon-oxygen core. This determines the 
rate of mass transfer to the secondary compan- 
ion and th e type of t h e rem nant left behind. We 
will follow IWebbinkl (|l979h in treating the first 
mass exchange modes for normal binary compo- 
nents, whose scheme accounts for the physical 
state of the star in more detail than the simple 
types of mass exchange (A, B C) introduced by 
Kippenhahn & Weigert 1 (|l967h . We will use both 
notations A, B, C and D (for very wide systems 
with independently evolving companions) for evo- 
lutionary types of binary as a whole, and Web- 
bink's notations for mass exchange modes for each 
component separately [la], [lb], [Ha], [lib], [Ilia], 
etc. 

3. Basic evolutionary states of normal 
stars 

The evolution of a binary system consisting ini- 
tially of two zero-age main sequence stars can be 
considered separately for each components until a 
more massive (primary) component fills its Roche 
lobe. Then the matter exchange between the stars 
begins. 

The evolutionary states of normal stars will be 
denoted by Roman figures (I-IV), whereas those 
of compact stars will be marked by capital letters 
(E, P, A, SA ...). We divide the evolution of a nor- 
mal star into four basic stages, which are signifi- 
cant for binary system evolution and bear a clear 
physical meaning. We will implicitly express the 
mass and radius of the star and the orbital semi- 
major axis in solar units (to = M/Mq, r = R/Rq, 
a = A/Rq), the time in million years, the lumi- 
nosities in units of 10 38 erg s _1 , the wind velocities 
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in units of 10 8 cm s _1 and the accretion rates M 
onto compact objects in units of lCP 8 A/0 yr~ 4 , 
unless other units are explicitely used. 

3.1. Main sequence stars 

At this stage, the star is on the zero-age main 
sequence (ZAMS) and its size is much smaller than 
the Roche lobe radius. The time the star spends 
on the main sequence is the core hydrogen burning 
time, tn~ which depend s on the stellar mass only 
(llben fe Tutukov Ill987l ): 



1.0 + 0.95^, 

103-9 — 3.8 log m+log m 

2400m" 



-2.16 



m > 79.0, 
79 > m > 10, 
10 > to > 2.3, 
to < 2.3, 



(6) 



10 4 m- 3 - 5 , 

The radius of the ZAMS star is assumed to be 



10 0.661ogm+0.05 ) m > ]_ 12 , 
TO, 771 < 1.2, 



(7) 



and its luminosity is 



log L = < 



-5.032 
-4.253 
-4.462 
-3.362 
-3.636 



2.65 log to, 

4.81ogTO, 

3.81ogTO, 

3.01ogTO, 

2.71ogTO, 



(a) 

(7) 
(<*) 
(*) 



(8) 



here we assume the next indication: (a), to < 0.6; 
(/3), 0.6 < m < 1.0; (7), 1.0 < to < 10.0; (S), 
10 < to < 48.0; (e), m > 48.0. 

The initial mass of the primary and the mode 
of the first mass exchange (which is determined by 
the initial o rbital period and masses of the com- 
ponents; see lWebbink l (fl979h determine the mass 
and the type of the core that will be formed during 
stage I. For example, for single stars and primaries 
in "type C" binaries that fill its Roche lobe having 
a degenerate core, we use the expressions 



m n = < 



{) .\jTl rflax , 

0.446 + 0.106m, 

min(0.36TO - 55 

0A4m maxJ 

~ M Ch , 
Im 14 



44to - 42 ) 

max, ^ "max / 



(a) 

(/J) 
(7) 
(S) 
(e) 
(0 
(V) 



(9) 



where m ma x is the maximum mass the star had 
during the preceding evolution. We assume the 
next indication in this formula: (a), m max < 0.8 
09), 0.8 < m max < 2.3; (7), 2.3 < m max < 4.0 
(6), 4.0 < m max < 7.5; (e), 7.5 < m max < 8.8 
(C), 8.8 < m rnax < 10.0; (77), m max > 10.0. 

A main-sequience star accreting matter during 
the first mass transfer will be treated as a rapidly 
rotating "Be-star" with the stellar wind rate dif- 
ferent from what is expected from a single star of 
the same mass (see below). 

3.2. Post main-sequence stars 

The star leaves the main sequence and goes to- 
ward the red (super)giant region. The star still 
does not fill its Roche lobe. The duration of this 
stage for a binary component is not any more a 
function of the stellar mass only (as in the case 
of single stars), but also depe nds on the initial 
binary type (A, B, or C) (see llben fc Tutukov 
(|l985Ul987h ;): 



tn = 



0, 


(«), 




(0), 


6300m~ 3 - 2 , 


(7), 


tHe, 


(*)> 



(10) 



In type A systems, the primaries fill their Roche 
lobes when they belong to the main sequence. We 
assume the next indication in this formula: (a), 
type A; (/3), type B excluding mode [IIIA]; (7), 
types C, D and mode [IIIA], m max < 5; (5), types 
C, D and mode [IIIA], m max > 5. 

The radius of the post-MS star rapidly increases 
(on the thermal time scale) and reaches the char- 
acteristic giant values. The star spends the most 
time of helium burning with such large radius. In 
the framework of our approximate description, we 
take the radius of the giant star to be equal to 
the maximum value, which depends strongly on 
the mass of its core and is calculated according 
to Webbink's mass tr a nsfer modes as follows (see 
also llben fc Tutukov I (|l985l . ll987h ): 



3000m 4 , 
1050(m c 
10m° 44 , 



0.5) 



0.(iS 



(a) 

(/?) 
(7) 



(11) 



We assume the next indication in this formula: 
(a), mode [IIIA] or [IIIB] with He core; (/3), 
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mode [IIIB] with CO or ONeMg core; (7), modes 
[I] or [II]. Formula (jTTJ) depicts stars with mass 
< 10M S . 

This maximum radius can formally exceed the 
Roche lobe size; in such cases we put it equal 
to 0.9i?L during the stage II. The most sensi- 
tive to this crude approximation are binaries with 
compact companions, which can lead, for exam- 
ple, to the overestimation of the number of ac- 
creting neutron stars observed as X-ray pulsars. 
However, these stages are less important for our 
analysis than the stages at which the optical star 
fills its Roche lobe. A more detailed treatment 
of normal star e v olutio n (given, for example, by 
Pols &; Marinus I (|1994l )) can reduce such uncer- 
tainties. 

Lum inosities of giants are taken from lde Jager I 



(1980): 



log/77 



15.92ro° 

(1. 0+m*)(2. 512+3. 162m e ) ' 
-^Q-4.462+3.81ogm 

1 q-3. 362+3.0 log m 

IQ-3. 636+2. 71ogrn 



(a) 

(7) 
(S) 



(12) 



We assume the next indication in this formula: 
(a), m < 23.7, m c < 0.7; (J3), to < 23.7, m c > 0.7; 
(7), 48 > m > 23.7; (5), to > 48. 

Radii of (super)giants are determined by us- 
ing the effective temperature T e ff and luminosi- 
ties. Typical effective temperatures are taken from 



Allen I (|l973h : 




(a) 

08) 
(7) 



(13) 



We assume the next indication in this formula: 
(a), m > 10.0; (/3), to < 10.0, type C or D; (7), 
to < 10.0, type A or B. So, we calculate .R77 with 
mass higher than IOMq calculate using formula 



r n = exp {2.3(0.5 log l u - 21ogT e// + 9.7)}. 



(14) 



3.3. Roche lobe overflow 



At this stage the star fills its Roche lobe (RL) 
and mass transfer onto the companion occurs. 
The mass transfer first proceeds on the thermal 



time scale (see extensive dis c ussion of this approx- 



imaiton in van den Heuvel (1994)) 



t KH ~ XknW/Le)- 1 , 



(15) 



The common envelope stage (CE) may be 
formed if the Roche lobe overflow occurs in the 
type C system (where the primary has a well de- 
veloped core) even for q < 1; otherwise (for type 
B systems) we use the condition q < q cr = 0.3 for 
the CE stage to occur. Radius of the star at the 
Roche lobe filling stage is tak en to be that of the 
equivalent Roche lobe radius (Eggleton I ll983h : 



Ri 



0A9q 



2/3 



o(l -e) 0.6g 2 / 3 + ln(l + q 1 / 3 )' 



(16) 



Here a is the binary orbital separation and e is the 
orbital eccentricity, q is the arbitrary mass ratio. 

For q < 0.6 a more precise approximation can 
be used: 



Rl f q 

— = 0.4622 — — 
a \l + q 



1/3 



(17) 



A star filling its Roche lobe has quite differ- 
ent boundary conditions in comparison with sin- 
gle stars. The stellar radius at this stage is limited 
by the Roche lobe. If the stellar size exceeds the 
Roche lobe, the star can lose matter on a time 
scale close to the dynamical one until its radius 
becomes smaller than the new Roche lobe size. 

Now we consider how the RL-filling star loses 
matter. Let the star be in equilibrium and 
R eq (M) = Rl(M) at the initial moment of time. 
When a fraction of mass 8m is transported to the 
companion, the mass ratio q and semi-major axis 
a of the binary changes depending on the mass 
transfer mode assumed (see below for details). 
The RL size then becomes equal to Rl(M — 5m). 

On the other hand, mass loss disturbs equi- 
librium of the star (hydrodynamical and ther- 
mal) . The hydrodynamical equilibrium is restored 

on the dynamical timescale td oc (GM/i? 3 ) 
The stellar radius changes to a value R a d(M — 
5m) (where "ad" means adiabatic), which can 
be bigger or smaller than the equilibrium radius 
R eq (M — 8m) of the star. The thermal equilib- 
rium establishes on the thermal time scale Tkh ~ 
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GM 2 /Rl, so after that the stellar radius relaxes 
to the equilibrium value R eq {M — 5m). 

Relations Rl(M — 5m), R a d{M — 5m) and 
R eq (M — Sm) determine the mode of ma ss transfer 
during the RL overflow stage. Following I Webbinkl 
( 19851 ), one usually introduces the logarithmic 
derivative C = dhxR/dhiM (R oc M c ). It lo- 
cally fits the real dependence R(M). Three values 
of Q are relevant: 



a 

Cad 
Ceq 



din R L 
dlnM ' 
d In R ad 
dlnM ' 
d In R eq 
d\nM ' 



(18) 



where "L", "ad" and "eg" correspond to the values 
of radii discussed above. 

Three possible cases are considered depending 
on Q: 

1- If Cad < Cl, the star cannot be inside its RL 
regardless of the mass loss rate (dM/dt < 0). 
Such stars lose their matter in hydrodynami- 
cal time scale. Mass loss rate is limited only 
by the speed of sound near the inner La- 
grangian point L%. Ceq is unimportant be- 
cause the size of the star becomes bigger and 
bigger than R^. The equilibrium is impossi- 
ble. 

2. Ceq < Cl < Cad- The star losing mass cannot 
be in thermal equilibrium, because otherwise 
its size would exceed Rl- Nevertheless, in 
this case R a d < Rl- So the hydrodynamical 
equilibrium is established. As a result, the 
star loses mass on thermal time scale. 

3. Cl < Cad, Ceq- In this case the size of the star 
losing mass becomes smaller than its RL. 
The evolutionary expansion of the star or 
the binary semi-major axis decrease due to 
orbital angular momentum loss via magnetic 
stellar wind (MSW) or gravitational radia- 
tion (GW) support the permanent contact of 
the star with RL. The star then loses mass 
on a time scale dictated by ots own evolu- 
tionary expansion or on a time scale corre- 
sponding to the orbital angular momentum 
loss. 



For non-degenerate stars, R eq increases mono- 
tonically with M. On the other hand, the expo- 
nent Cad is determined by the entropy distribution 
over the stellar radius which is different for stars 
with radiative and convective envelopes. It can be 
shown that stars with radiative envelopes should 
shrink in response to mass loss, while those with 
convective envelopes should expand 5 . 

Therefore, stars with convective envelopes in bi- 
naries should generally have a higher mass loss 
rate than those with radiative envelopes under 
other equal conditions. The next important fac- 
tor is the dependence Rl{M). It can found by 
substituting one of the relations a(M) (see below) 
into equation (|16p or into equation (117[) and dif- 
ferentiating it with respect to M. For example, 
assuming the conservative mass exchange when 
the total mass and the orbital angular momen- 
tum of the system do not change, one readily gets 
that the binary semi-major axis decreases when 
the mass transfer occurs from the more massive to 
the less massive component; Rl decreases of the 
primary correspondingly. When the binary mass 
ratio reaches unity, the semi-major axis takes on a 
minimal value. In contrast, if the less massive star 
loses its mass conservatively the system expands. 
In that case the mass transfer can be stable. 

If more massive component with radiative en- 
velope fills its RL, the mass transfer proceeds on 
thermal time scale until the masses of the compo- 
nents become equal 6 . The next stage of the first 
mass exchange poceeds in more slower (nuclear) 
time scale. If the primary has convective enve- 
lope, the mass transfer can proceed much faster 
on a time scale intermediate between the thermal 
and hydrodynamical one, and probably on the hy- 
drodynamical time scale. In that case the fast 
stage of the mass exchange ends when the mass 
of the donor decreases to ~ 0.6 mass of t he sec - 

1978h . 



ondary companion (jTutukov fc Yungelson 



Further mass transfer should proceed on the evo- 
lutionary time scale. 



5 The adiabatic convection in stellar envelope can be de- 
scribed by the polytropic equation of state P oc p 5 / 3 , simi- 
lar to non-relativistically degenerate white dwarfs. For such 
equation of state the mass-radius relation becomes inverse: 
R oc M -1 / 3 . For non-degenerate stars with convective en- 
velopes this relation holds approximately. 

6 The mass exchange can stop earlier if the entire envelope 
is lost and the stellar core is stripped (the core has other 
values C ad and Ceq)- 
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The process of mass exchange strongly depends 
on stellar structure at the moment of the RL over- 
flow. The structure of the star in turn depends 
on its age and the initial mass. The moment of 
the RL overflow is determined by the mass of the 
components and by the initial semi-major axis of 
the system, to calculate a diagram in the M — a 
(or M — P rb) plane which allows us to conclude 
when the primary in a binary with given initial pa- 
rameters fills its RL, what is its structure at that 
moment and what type of the first mass exchange 
is expected. W e use the diagram calculated by 
Webbink ( 1979h (see also the description of mod- 
ern stellar wind scenarios below). 

We distinguish different sub-stages of the RL 
overflow according to the characteristic timescales 
of the mass transfer: 

stage III: 

This is the most frequent case for the first mass 
transfer phase. The primary fills its RL and the 
mass transfer proceeds faster than evolutionary 
time scale (if outer layers of the star are radia- 
tive, then it is thermal time scale, if outer layers 
are convective, then time scale is shorter, up to 
hydrodynamical time scale). This stage comes to 
the end when mass ratio in the system changes 
("role-to-role transition"), i.e. when the mass of 
the donor (mass losing) star is equal to mass of 
second companion (for radiative envelopes) or 0.6 
of the mass of the second companion (for convec- 
tive envelopes). This stage also stops if the donor 
star totally lost its envelope. 

stage Hie: 

This is the slow (evolutionary driven) phase of 
mass transfer. We assume it to occur in short- 
period binaries of type A. However, it is not ex- 
cluded that it may occur after the mass reversal 
during the firs t stage of mass ex c hange for binaries 
of type B (see van den Heuvel ( 199J)), e.g. as in 
wide low-mass X-ray binaries. 

stage Ills 

This is the specific to super-accreting com- 
pact companions substage of fast mass transfer at 
which matter escapes from the secondary compan- 
ion carrying away its orbital angular momentum. 
Its duration is equal to 



tills = tKH' 



g(l + g) 
\-q-2q 2 



(19) 



q = M a /M d < 0.5, 

here and below subscripts "a" and "d" refer to 
the accreting and donating mass star, respectively. 
For systems with small mass ratios, q < 0.5, this 
timescale corresponds to an effective q-time short- 
ening of the thermal time for the RL-overflowing 
star. 

stages IIIm,g 

At these stages, the mass transfer is controlled 
by additional losses of orbital angular momentum 
J or b caused by magnetic stellar wind or gravita- 
tional wave emission. The characteristic time of 
the evolution is defined as t j = —(J or bl Jorb) , 
and in the case of MSW is (see Verbunt fc Zwaan 



(119811 ): llben fc Tutukov I (119871) 1 



T MSW 



4.42- 



(mi 



W2) 2 m 4 



(20) 



op 



Here m op denotes mass of the low-mass optical 
star (0.3 < m < 1.5) that is capable of producing 
an effective magnetic stellar wind (because only 
such stars have outer convective envelopes which 
are prerequisit for effecftive MSW), A is a numer- 
ical parameter of order of unity. We have used 
the mass-radius relation r rs m for main sequence 
stars in deriving this formula. The upper limit 
of the mass interval and empirical br aking law for 
main- sequence G-stars are taken from ISkumanich" 1 



(|1972T ). the lower limit is determined by absence of 



cataclysmic variables with 


orbital periods j 


°orb ~ 


3^ dVerbunt 19841: iMestel 


1952; 


Kawalcr 


19881: 


Tout & Prinde Il992l IZanerilli et al. Ill997 


). We 


use A = 1 (see for details Kaloeera & Webbink 



The time scale of the gravitational wave emis- 
sion is 



tgw = 124.2 



a 



(21) 



+ 24 + 96 ' 1 ' ' 



Wether the evolution is governed by MSW or 
GW is decided by which time scale {tmsw or 
tgw) turns out to be the shortest among all ap- 
propriate evolutionary time scales. 

stage Illwd 

This is a special case where the white dwarf 
overflows its RL. This stage is encountered for very 
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short period binaries (like Am CVn stars and low- 
mass X-ray binaries like 4U 1820-30) whose evolu- 
tion is controlled by GW or MSW. The mass trans- 
fer is calculated using the appropriate time scale 
(GW or MSW). The radius of the white dwarf in- 
creases with mass Rwd oc A/ -1 / 3 . This fact, how- 
ever, does not automatically imply that the mass 
transfer is unstable, since the less massive WD 
fills its RL first. It can be shown that the mass 
exchange is always stable in such systems if the 
mass ratio q < 0.8. This condition always holds 
in WD+NS and WD+BH systems. WD loses its 
matter until its mass decreases to that of a huge 
Jupiter-like planet (~ a few 1O _3 M0), where the 
COulomb interaction reverses the mass-radius re- 
lation R(M). Such a planet can approach the sec- 
ondary companion of the system due to GW emis- 
sion until the tidal forces destroy it completely. 
The matter of the planet can fall onto the sur- 
face of the second companion ir form a long-living 
disk around it. If the second star is a neutron star 
and its rotation had been spun up by accretion 
such that a millisecond radio pulsar appeared, the 
planet can be evaporat ed by relativistic particles 
emm i ted by the pulsar fPaczvnski fc Sienkiewicz 



Kalogera fc Webbink 



19981) . 



1983t IJoss fc Rappaport Ill983l : iKolb et al. Ill998: 



The mass loss rate at each of the Ill-stages is 
calculated according to the relation 



M = AM I n 



(22) 



where AM is the a priori known mass to be lost 
during the mass exchange phase (e.g. AM = Mi — 
(Mi -I- M<z)/2 in the case of the conservative stage 
III, or AM = Mi-M core (Mi) in case of III(e,m,g) 
or CE) and r, is the appropriate time scale. 

The radius of the star at stage III is assumed 
to follow the RL radius: 



Rf(M d (t)) = R d (M d (t)). 
3.4. Wolf-Rayet and helium stars 



(23) 



In the process of mass exchange the hydrogen 
envelope of the star can be lost almost completely, 
so a hot white dwarf (for m < 2.5), or a non- 
degenerate helium star (for higher masses) is left 
as a remnant. The life-time of the helium star is 
deter mined by the helium b urning in the stellar 
core (|lben fc Tutukov II 19851 ) 



tHe 



1658m" 2 , 


(a) 


1233m- 38 , 


id) 


0.1t H , 


(7) 


6913m- 347 , 


(S) 


-10, 


(e) 


O.U H , 


(0 



(24) 



We assume the next indication in this formula: 
(a), m < 1.1, modes [IIA-IIF]; (/?), m > 1.1, 
m max < 10, modes [IIA-IIF]; (7), m max > 10, 
modes [IIA-IIF]; (<5), mode [IIIA]; (e), mode [IIIC]; 
(C), modes [IIIB,D,E] and type D. 

If the helium (WR) star fills its Roche lobe 
(a r elatively rare so-called "BB " case of evolu- 
tion: iDelgado fc Thomas 1 ( 1981 ): see discussion in 
Ivan den Heuvel I (|l994h ). the envelope is lost and 
a CO stellar core is left with mass 



1.3 + 0.65(m - 2.4), m > 2.5, 
0.83m - 36 , m < 2.5, 



(25) 



The mass-radius depende nce in this case is 
(|Tutukov fc Yungelson ll 19731) 



o.(> 



r WR = 0.2m 



3.5. Stellar winds from normal stars 



(26) 



The effect of the normal star on the compact 
magnetized component is largely determined by 
the rate M and the velocity of stellar wind at in- 
finity Vqo, which is assumed to be 



= 3w p « 1.85\/ m/r, 



(27) 



where v p is the escape velocity at the stellar sur- 
face. 

For "Be-stars" (i.e. those stars at the stage "I" 
that increased its mass during the first mass ex- 
change), the wind velocity at the infinity is taken 
to be equal to the Keplerian velocity at the stellar 
surface: 



= sfGMjR S3 0.44 vmA 



(28) 



The lower stellar wind velocity leads to an effec- 
tive increase of the captured mass rate by the sec- 
ondary companion to such "Be-stars" . 

The stellar wind mass loss rate at the stage "I" 
is calculated as 
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(jde Jaeer ||l980l) 

m = 52.3a„,Z/u 0O , 



(29) 



Here a m = 0.1 is a numerical coefficient (in gen- 
eral, we can treat it as free parameter). 

For giant post-MS stars (stage "II" ) we assume 
Voa = 3w p and for massive star we take the max- 
imum between the stellar wind r ate given by de 
Jager's formula and that given bv iLamers" (1981) 



m = max(52.3a IU , 10 



2.33 ' 



Voo III 

M > 1OM ; 



0.0!) 



■). (30) 



For red super-giants we u se Reimers's formula 
(|Kudritzki fc Refiners |[l978h : 



I Ir 

rh = max(52.3a UI ,1.0 — ), (31) 

Woo m 

M > 10M Q ; 

For a Wolf-Rayet star the stellar wind loss rate 
can significantly increase (up to 1O _5 M year -1 ). 
We parametrize it as 



MwR — kwuMwil/tHei 



(32) 



where the numerical coefficient is taken to be 
kwR = 0.3 (in general, it can be changed if nec- 
essary). The mass loss in other stages (MS, (su- 
per )giant) is assumed to be limited by 10% of the 
mass of the star at the beginning of the stage. 

3.6. Change of binary parameters: mass, 
semi-major axis and eccentricity 

The duration of any evolutionary stage is de- 
termined by the more rapidly evolving compo- 
nent At — min(Aii, At^)- On the other hand, 
based on the evolutionary considerations we are 
able to calculate how the mass of the faster evolv- 
ing star changes (e.g. due to the stellar wind or 
RL overflow) , that is we can estimate the quantity 
AM = Mi — Mf. Then we set the characteristic 
mass loss rate at this stage as 



Mo = AM/ At, 



(33) 



Next, we should calculate the change of mass for 
the slower evolving companion, the orbital semi- 
major axis and the eccentricity. 



3.7. Mass change 

The mass of the star loosing matter is calcu- 
lated as 



Mf — Mi — M x At, (34) 
Accordingly, the mass of the accreting star is 

M f = Mi + M c x At, (35) 

where M c is the accretion rate of the captured 
matter. 

For stages without RL overflow the accretion 
rate of the captured stellar wind matter is 



m c = 3.8x10" 



a(vl +0.19(toi +m 2 )/a) 



(36) 

At stages where RL overflow ocurs and both 
components are normal (non-degenerate), we will 
assume that the accretor can accomodate mass at 
the rate determined by its thermal time, i.e. 

. / t KH (donor) \ 

m c = m [ — — . (37) 

\tKH [accretor) J 

This means that the evolution can not be 
fully conservative, especially during the first mass 
transfer where the primary component usually has 
a shorter thermal time scale. 

The mass increase rate by compact accretors 
is assumed to be limited by the critical Edding- 
ton luminosity (see, however, the possible hyper 
accretion stage discussed below) : 



AnGMrrir, , H , 

L E dd= ^« 1.3-10 38 xm erg/s (38) 

(<tt is the Thomson cross-section) at the stopping 
radius R s top for the accreted matter (see, e.g., de- 
tailed discussion in iLipunov I dl992h l. This corre- 
sponds to the critical accretion rate 



• _ R L Edd 

lVl cr — J^stop 



(39) 



Thus, the mass of the accreting compact star at 
the end of the stage is determined by the relation 

M f = Mi + min(M c , M cr ) x At . (40) 



9 



3.8. Semi-major axis change 

The binary separation a changes differently for 
various mass exchange modes. First, we introduce 
a measure of non-conservativeness of the mass ex- 
change as the ratio between the mass change of 
the accretor and the donor: 



P = -(Mi ccr -Mf ccr )/(Mi 



Minor)- (41) 



If the mass exchange is conservative (f3 = 1, 
i.e. M a + Md = const) and one can neglect the 
angular momenta of the components, the orbital 
momentum conservation implies 



M s a M{ 



(42) 



In a more general case of quasi-conservative 
mass transfer < j3 < 1, the orbital separation 
changes differently depending on the specific an- 
gular momentum carried away from the sys t em b 
the escaping matter (see Ivan den Heuvel I <|199 
for more detail). We treat the quasi-conservative 
mass transfer by assuming the isotropic mass loss 
mode in which the matter carries away the spe- 
cific orbital angular momentum of the accreting 
component (j a ) 



Jout = (1 - P)M c j a . 
From here we straightforwardly find 




3+2//3 



(43) 



(44) 



In this case the orbital separation always increases. 

When the orbital angular momentum is carried 
away by GW or MSW with no RL overflow, the 
following approximate formulas are used: 



(1- M/tmswY^M MSW, 
(1- M/rawf^M GW, 



(46) 



In a special case of a white dwarf filling its 
RL (stage "Illwd" above), assuming a stable (i.e. 
where dh\R w d/ d\xiM = dh\ Rb,l I d\n. M) conser- 
vative mass transfer with account for the mass- 

— 1/3 

radius relation R w d oc M wd ' , the orbital separa- 
tion must increase according to the equation 



C 1L 

a, 



mj_ 

rrii 



-2/3 



(47) 



where rrn is the initial mass of the WD donor and 
rrif is its mass at the end of the mass transfer. 

3.9. The change of eccentricity 

Tidal interaction between components, as well 
as the orbital angular momentum loss due to MSW 
or GR decrease the eccentricity of the binary sys- 
tem. The tidal interaction is essential in very 
close binaries or even during the common envelope 
stage. MSW is effective only in systems with low- 
mass late- type main sequence stars (see above), 
GW losses become significant only in short-period 
binaries. 

The tidal interaction conserves the orbital an- 
gular momentum which implies the relation 



In this formula q — M accr /Md on or and the 
non-conservative parameter (3 is set to be the 
minimal value between f3 = 1 and the ra- 
tio Tk h {donor) /TKH{accr) {Tkh {donor) and 
TKn{accr) is the thermal time of donor and ac- 
cretor, respectively). 

When no matter is captured by the secondary 
companion without additional losses of angu- 
lar momentum (the so-called "absolutely non- 
conservative case"), which relates to the spherical- 
symmetric stellar wind from one component, we 
use another well-known formula 



a f _ M\ + M l 2 



Mi + M, 



f 



(45) 



(1 — e 2 ) = const. 



(48) 



It seems that the orbit becomes a circle faster 
than major semi-axis of the orbit decreases during 
common envelope stage. We suppose that t cyr — 
l/3tcB- We accept that t cyr for RL-filling stars 
is equal to its Kelvin-Hclmholtz time. Detached 
systems change their eccentr i city during the next 
character time ( Zah nl Il975t iPress fc T eukolskv 
1977t IZahn Ill989bt IZahn fc Bouchet Ill989bf ) 



tKH 



1 - e 



3/2 



R_ 



(49) 
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here tKH ~ GM 2 /Rl is thermal time of the star, 
R is radius of the star, Rl is its RL size. For sys- 
tems which consist of two normal stars we choose 
minimal value of t cyr . Resonances at ve ry high ec 



centric ities are not taken into account (jMardling 
ll995alTbh . 

Orbits of the systems with MSW become circu- 
lar during t cyr = t msw (see HU) . 

In case of GW analytical exact solutions were 
obtained for a(t) and e(i) (jPeters k, Mathews 



19631 : IPeters lll964h 



4. Special cases: supernova explosion and 
common envelope 

Supernova explosion in a binary is treated as 
an instantaneous mass loss of the exploding star. 
The additional kick velocity can be imparted to 
the newborn neutron star due the collapse asym- 
metry (see below for discussion). In this case the 
eccentricity and semi-major axis of the binary af- 
ter th e explosion can be straightforwardly calcu- 
lated (iBoersma lll96ll) (see necessary formulas also 
in Grishchuk et al. 2001). Briefly, we use the fol- 
lowing scheme. 

1. First, velocities and locations of the compo- 
nents on the orbit prior to the explosion are 
calculated; 

2. then the mass of the exploding star M pr - 
M remnan t is changed and the arbitrarily di- 
rected kick velocity w is added to its orbital 
velocity; 

3. after that the transition to the new system's 
barycenter is performed (at this point the 
spatial velocity of the new center of mass of 
the binary is calculated); 

4. in this new reference frame the new total en- 
ergy E' tot and the orbital angular momentum 
J' orb are computed; if the new total energy is 
negative, the new semi-major axis a' and ec- 
centricity e' are calculated by using the new 
J'orb en d E'tot i ^ the total energy is positive 
(that is, the binary is unbound) spatial ve- 
locities of each component are calculated. 



The kick velocity w distribution is taken in the 
Maxwellian form: 



ft \ w 

f(W) 2 e 



(50) 



We suppose that the absolute value of the ve- 
locity that can be added during the formation of 
a black hole depends on the mass loss by the col- 
lapsing star, the value of the parameter wq during 
the BH formation is defined as 



w b Q h = (1 - k bh )w . 



(51) 



An effective spiral-in of the binary compo- 
nents occurs during the common envelope (CE) 
stage. This com plicated process (introduced by 



Paczvnski I ([1976)) is not fully understood as yet, 



so we use the conventional energy consideration 
to find the binary system parameters after the CE 
by introducing a parameter olce = AE b j AE orb , 



where AEh 



E n 



Ethe 



is the binding en- 



ergy of the ejected envelope matter and AE or i, is 
the drop in the orbit al energy of the system during 



the spiral-in phase (jvan den Heuvel 1 119941 ). This 



parameter measures the fraction of the system's 
orbital energy that comes during the spiral-in pro- 
cess to the binding energy (gravitational minus 
thermal) of the ejected common envelope. Thus 



acE 



( GM a M c GM a M d \ GM d (M d -M c ) 



V 2a, 



2a< 



R d 



(52) 

where M c is the mass of the core of the mass- 
losing star with the initial mass M d and radius R d 
(which is simply a function of the initial separation 
di and the initial mass ratio M a /M d , where M a is 
the mass of the accreting star). 

On the CE stage the luminosity of the accret- 
ing star can reach the Eddington limit so that the 
further increase of the accretion rate can be pre- 
vented by radiation pressure. This usually hap- 
pens at a ccretion rat e s M ~ 10~ 4 — 1O _5 M0 yr" 1 . 
However, IChevalier (1993) suggested that when 
the accretion rate is higher (M ~ 10~ 2 — 10 -3 M Q 
yr -1 ), the energy is radiated away not by high- 
ene rgy photons only , but also by neutrinos (see 
also Zeldovich et al. ( 19721 ) and the next section). 
On the typical time scale for the hyper accretion 
stage of 10 2 yr, up to ~ IMq of matter may be 
incident onto the surface of the neutron star. 
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5. Three regimes of mass accretion by 
neutron stars 

A considerable fraction of observed neutron 
stars have increased their masses in the course of 
their evolution, or are still increasing their masses 
(e.g., in X-ray sources). But how large can this 
mass increase be? It is clear that the only origin 
of a mass increase is accretion. It is evident that 
the overall change in the mass of a neutron star 
is determined not only by the accretion rate, but 
also by the duration of the accretion stage: 



AM 



Mdt = MT a , 



(53) 



where M is the mean accretion rate and T a is the 
lifetime of the accretion stage. We emphasize that, 
in the case under consideration, the accretion rate 
is the amount of matter falling onto the surface 
of the neutron star per unit time, and can differ 
signifcantly from the values indicated by the clas- 
sical Bondi-Hoyle formulas. Three regimes of ac- 
cretion are possible in a close binary containing a 
neutron star: ordinary accretion, super-accretion, 
and hyper-accretion. 

5.1. Ordinary accretion 

The ordinary accretion regime is realized when 
all matter captured by the gravitational field of the 
neutron star falls onto its surface. This is possible 
only if the radiation pressure and electromagnetic 
forces associated with the magnetic field of the star 
and its rotation are small compared to the gravita- 
tional force. In this case, the increase in the mass 
will be precisely determined by the gas dynamics 
of the accretion at the gravitational-capture ra- 
dius or, if the donor fills its Roche lobe, by the 
binary mass ratio and the evolutionary status of 
the optical component. In this case, the accretor 
is observed as an X-ray source with luminosity 



L x = M 



GM X 



(54) 



where M x and i?* are the mass and radius of the 
neutron star. The accretion rate M is determined 
by the Bondi-Hoyle formula 



where Rq is the gravitational-capture radius of 
the neutron star, v is the velocity of the gas flow 
relative to the neutron star, and p is the density 
of the gas. 

The X-ray luminosity of the accretor L x and 
its other main parameters can be used to estimate 
the mass AM accumulated during the accretion 
phase: 

A " - (56) 

5.2. Super-accretion 

Regime of supe r accretion was c onsidered, for 
instance, in paper iLipunov I (|1982d). Despite the 
absence of detailed models for supercritical disk 
accretion (supercritical accretion is realized pre- 
cisely via an accretion disk), it is possible to esti- 
mate the main characteristics of the process - the 
accretion rate, magnetosphere radius, and evolu- 
tion equations. Accretion is considered to be su- 
percritical when the energy released at the radius 
where the accretion exceeds the Eddington limit: 



M— — - > L Edd 



R 



stop 



1.ZS>x\Q m {M x /Mq) ergs" 1 , 

(57) 

where R s top is either the radius of the neutron star 
or the magnetosphere radius Ra- 

For strongly magnetized neutron stars with 
magnetic fields B 3> 10 s G, all matter arriving 
at the magnetosphere is accreted onto the mag- 
netic poles, where the corresponding gravitational 
energy is released. If the black body temperature 
T, roughly estimated as 



SaT 4 = M 



GM X 
1 I 



(58) 



is higher than 5 x 10 9 K (S is the area of the 
base of the accretion column) , most of the energy 
will escape from the neutron star in the form of 
neutrinos , and, hence, will not hinder accretion 



(|Zeldovich et al. Ill972t iBasko fc Sunvaev Ill975h . 
In this case, the rate at which the neutron star 
accumulates mass will be 



M ~ M r . 



M = irR G pv . 



(55) 



Ra 
R* 



> M c . 



(59) 
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For lower temperatures there should be an up- 
per limit on the accretion rate equal to the stan- 
dard Eddington limit. 

5.3. Hyper-accretion 

A considerable fraction of neutron stars in bi- 
nary systems pass through the common-envelope 
stage in the course of their evolution. In this case, 
the neutron star is effectively immersed in its opti- 
cal companion, and for a short time (10 2 -10 4 yrs) 
spirals-in inside a dense envelope of the compan- 
ion. The formal accretion rate estimated using 
the Bondi-Hoyle formulas is four to six orders of 
magnitude higher than the critica l rate and, as 
was suggested by IChevalier ( 1993 ). this may re- 
sult in hyper-accretion, when all the energy is car- 
ried away by neutrinos for the reasons described 
above. There are currently no detailed theories 
for hyper-accretion or the common-envelope stage. 
The amount of matter accreted by the neutron 
star can be estimated as 



AM = l ThyPerl -^] 2 Mdt 



(60) 



1 



(M opt - M core ) 



( M x 



Op I 



where Th VP er is the duration of the hyper- accretion 
stage, Rg is the gravitational-capture radius of 
the neutron star, a is the initial semi-major axis 
of the close binary orbit, M core is the mass of the 
core of the optical star, and M opt and M x are the 
total masses of the optical star and of the neu- 
tron star at the onset of the hyper- accretion stage. 
The mass of the neutron star can increase during 
the common envelope st age as much as ~ IMq 



(jBogomazov et al. 1 120051) . up to Mqv- Such NSs 



collapse into black holes. 

6. Mass accretion by black holes 

If the black hole has formed in the binary sys- 
tem, its X-ray luminosity is 



L, 



(61) 



where /i = 0.06 for Schwarzschild black hole and 
/i = 0.42 (maximum) for extremal Kerr BH. 

We use the Bondi-Hoyle formulas to estimate 
the accretion rate onto BH. 



Powerful X-ray radiation is able to originate 
only if an accretion disc has formed around the 
black hole (jKarpov fe Lipunovl 120011 ). For spher- 
ically symmetric accretion onto a black hole the 
X-ray luminosity is insignificant. A very low stel- 
lar w ind velocity is n ecessary to form an accretion 
disc (jLipunov Ill992h 



V < V cr « (62) 

« 320(4r / ) 1 / 4 m 3 / 8 T 10 1/4 ^ 1/8 (l + tan 2 /3)- 1 / 2 , 

where r\ is averaged over the z-coordinate dynamic 
viscous coefficient, m = M x /Mq, M x is the rela- 
tivistic star mass, Tio = T/10, T is the orbital 
period in days, Rg — i? m i„/10 8 cm, R m i n is the 
minimal distance from the compact object up to 
which free Keplerian motion is still possible and 
(3 is the accretion disk axis inclination angle with 
respect to the radial direction. For black holes 
R mm = 3i? 9 , where R g = 2GM bh /c 2 . 

7. Accretion induced collapse and com- 
pact objects merging 

WD and NS are degenerate configurations, 
which have upper limit of their mass (the Chan- 
drasekhar and Oppenheimer-Volkov limits corre- 
spondingly). The Chandrasekhar limit depends 
on chemical composition of the white dwarf 



1.44M©, 
!/,•„ = { 1.40M©, 
1.38M , 



He WD, 
CO WD, 
ONeMg WD, 



(63) 



If the mass of WD becomes equal to Mch, 
the WD loses stability and collapses. The col- 
lapse is accompanied by the powerful thermonu- 
clear burst observed as a type la supernova. Col- 
lapses of He and CO W Ds leave no remnants 



(jNomoto k, Kondo II 199 11 ). The outcome of the 



collapse of a ONeMg WD is not clear. It can lead 
to the formation of a neutron star (the accretion 
induced collapse, AIC). Some papers come to the 
differ ent conclusion about t he res u lt of AICs (see 
e.g. iGarcia-Berro fc Iben ( 1994f ): Ritossa et al. 
(<1996l) 1. The question about the NS formation 



during the WD collapse remains undecided. Nev- 
ertheless, some NSs could hav e been formed from 
AIC WD (|van Paradijs Ill997l ). In the "Scenario 
Machine" code the possibility of NS formation 
during ONeMg WD is optional. 
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The merging of two compact objects in a binary 
WD system (e.g., due to the GW losses) should 
likely to be similar to AIC. During the merging 
of two helium WDs, one object with a mass of 
less than Mch (for example, 0.5M© + O.6M -> 
1.1 Mq) can form. At the same time, if the total 
mass of the components exceeds 1.0 — 1.2M , a 
thermonuclear burning can happen. It is likely 
that the merging of a ONeMg WD with another 
WD can form a NS. 

For the typical NS mass « 1.4M , the binary 
NS+NS merging event can produce a black hole. 
Massive (ss 2.8M ) neutron star can be formed 
only if the NS equation of state is very hard and 
M v - 2.8 - 3.OM . 

BH+BH merging should produce a rapidly ro- 
tating black hole with the mass equal to the total 
mass of the coalescing binary. 

8. Additional scenarios of stellar wind 
from massive stars 

8.1. Evolutionary scenario B 

In the end of 1980s and in the beginning of 
1990s the series of new evolutionary tracks were 
calc ulated. The authors used new tables of opaci- 
ties (jRogers fc Iglesias Ill991[ iKurucz Ill99ll). new 



cross-sections in nuclear reactions (Landre et al 



1990 ) and new paramet ers of convection in stars 
(|Stothers fc Chin I [l99lh . For stars with M < 
IOMq those tracks proved to be almost coincident 
with previous calculations. More massive stars 
had much stronger stellar winds. 

A massive star loses up to 90% of its ini- 
tial mass in the main-sequence, supergiant, and 
Wolf-Rayet stages via stellar wind. There- 
fore, the presupernova mass in this case can 
be 8 — IOMq, essential ly independent of the 



initial mass of the star rtde Jaeer et al.~ 1988t 



Nic uwenhuiizen &; de Jager 1 1 10901 : ISchaller et al. 
19921) . 



8.2. Evolutionary scenario C 

The papers mentioned above were criticised and 
in 1998 a new version of the evolutiona ry scenario 
was developed (|Vanbevern et al. 1119981 ) . The stel- 



lar wind loss rates were corrected taking into ac- 
count empirical data about OB and WR stars. 
Here we list the main equations of this scenario 



(all results concern only with the stars with initial 
mass M > 15M ). 



logM 



1.671ogL- 1.55 log T eJf -6 
log L + log R - log M - 7.5, 
0.81ogL-8.7, (7) 
log L- 10, (6) 



.29, (a) 

03) 

(64) 

We assume the next indication in this formula: 
(a), H burning in the core; (/3), giant, Mo > 
4OM ; (7), giant, M < 4OM ; (6), Wolf-Rayet 
star. 

Note that in this subsection we used the mass 
M is in Mq , the luminosity L is in L and the ra- 
dius R is in Rq. With these new calculations, the 
Webbi nk diagram described abo ve changed signif- 
~ (|l998h . 



icantly IVanbevern et al. 



In this scenario, the total mass loss by a star is 
calculated using the formula 



AM = (M — M core ), 



(65) 



where M core is the stellar core mass (l66l). If the 



maximum mass of the star (usually it is the initial 
mass of a star, but the mass transfer in binary 
systems is able to increase its mass above its initial 
value) M max > 15M , the mass of the core in the 
main sequence stage is determined using (|66b ). 
and in giant and supergiant stages using (1661 3). In 
the Wolf-Rayet star stage (helium star), if Mwr < 
2.5Mq and M max < 2OM it is described using 
(j66V ). if M WR > 2.5M and M max < 2OM as 
@5P), if M max > 2OM as ((66k) 



1.62m° p 8 t 3 

2Q-3.051+4.211gm o 

0.83m^ 3 | 
1.3 + 0.65(mwR - 
mcore — 3.03m' 



0.342 
opt 



(a) 

-0.93(lgm opt ) 2 
(7) 

2.4) (5) 



(e) 



(66) 

These evolutionary scenarios have some pecu- 
liar properties with respect to the classical sce- 
nario. One of them is that the strong stellar wind 
from massive stars leads to a rapid and significant 
increase of the system's orbit and such stars can- 
not fill its RL at all. 

There are three observational facts that conflict 
with the strong stellar wind scenarios: 
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1. A very high M is a problem by itself. The 
observers calculate this quantity for most 
of OB and WR stars using the emission 
measure EM cx J n 2 e dl. The estimate of 
M using EM is maximal for homogenous 
wind. However, there are e videnc es (see 



e.g. Cherepashchuk et al. (1984)) that 



stellar winds of massive stars are strongly 
"clumpy". In this case the real M must 
be 3-5 times less. This note is especially 
important for the scenario with high stellar 
wind. 

Very massive Wolf-Rayet stars do exist. 
There are at least three double WR+OB sys- 
tems including very massive WR-stars: CQ 
Cep 40 M^, HP 311884 48M W and HD92740 



77M Q (|Cherepashchuk et al. I fl996h . Such 



heavy WR stars are at odds with the as- 
sumed high mass loss rate. 

In the semi-detached binary system RY Set 
(W Ser type) the ma ss of the primary com 
pone nt is w 35M Q (jCherepashchuk et al. 
1996^. This mass is near the limit 



(|Vanbevern et al~ 19981 ) beyond which the 
star, according to the high mass loss sce- 
nario, cannot fill its RL. 

8.3. Evolutionary Scenario W 

The evolutionary scenario W is based on the 
stellar evolution calcu lations by Woosley et al. 
(jWooslev et al. I (|2002h . Fig. 16), which repre- 
sents the relationship between the mass of the 
relativistic remnant and the initial mass of the 
star. We included into population-synthesis code 
two models with W-type stellar winds, which we 
label Wb and Wc. In models Wb and Wc, the 
mass-loss rates were computed as in scenario B 
and scenario C, respectively. The use of these 
models to calculate the wind rat e in a scenario 
based on Woosley's diagram (jWooslev et al. 
(|2002h . Fig. 16) is justified by the fact that 



scenarios B and C are based on the same nu- 
merical express i ons fo r the mass-loss rat e s from 



Sch aller et al~l (fl9 92) 



Nieuwenhuiizen fc de J agcr 



used by I Woosley et al. I (|2002h . 



Vanbevern et al. 
(|l990h ~ 



(11998ft : 
that were 



9. The "Ecology" of Magnetic Rotators 

One of the most important achievements in as- 
trophysics in the end of the 1960s was the real- 
ization that in addition to "ordinary" stars, which 
draw energy from nuclear reactions, there are ob- 
jects in the Universe whose radiation is caused by a 
strong gravitational and magnetic field. The well- 
known examples include neutron stars and white 
dwarfs. The property that these objects have in 
common is that their astrophysical manifestations 
are primarily determined by interaction with the 
surrounding matter. 

In the early 1980s, this approach led to the cre- 
ation of a complete classification scheme involv- 
ing various regimes of interaction between neu- 
tron stars and their environment, as well as to the 
first Monte Carlo simulation of the NS evolution 



(jLipunov I Il984h . In addition to NSs, this scheme 
has been shown to be applicable to other types of 
magnetized rotating stars. 

By virtue of the relationship between the grav- 
itational and electromagnetic forces, the NS in 
various states can manifest itself quite differently 
from the astronomical point of view. Accordingly, 
this leads to the corresponding classification of NS 
types and to the idea of NS evolution as a gradual 
changing of regimes of interaction with the envi- 
ronment. The nature of the NS itself turns out to 
be important also when constructing the classifica- 
tion scheme. This indicates that there should be a 
whole class of quite different objects which have an 
identical physical nature. To develop the theory 
describing properties of such objects (in a sense, it 
should establish "ecological" links between differ- 
ent objects), it proved to be convenient to use sym- 
bolic notations elaborated for the particular case 
of NS. We start this subsection with recollecting 
the magnetic rotator formali sm (mainly according 
to the paper Lipunov ( 1987h ). 



9.1. A Gravimagnetic Rotator 

We call any gravitationally-bounded object 
having an angular momentum and intrinsic mag- 
netic field by the term "gravimagnetic rotator" 
or simply, rotator. In order to specify the intrin- 
sic properties of the rotator, three parameters are 
sufficient - the mass M, the total angular mo- 
mentum J — Iui (I is the moment of inertia and u> 
is the angular velocity), and the magnetic dipole 
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moment fi. Given the rotator radius Rq, one can 
express the magnetic field strength at the poles 
Bo by using the dipole moment Bo — The 
angle between the angular moment J and the 
magnetic dipole moment fi can also be of impor- 
tance: (3 = arccos(J/x). 

9.2. The Environment of the Rotator 

We assume that the rotator is surrounded by 
an ideally conductive plasma with a density 
and a sound velocity Ooo at a sufficiently far dis- 
tance from the rotator. The rotator moves rela- 
tive to the environment with a velocity Woo. Un- 
der the action of gravitational attraction, the sur- 
rounding matter should fall onto the rotator. A 
rotator without a magnetic field would capture a 
stationary flow of matter, M c , which can be es- 
tim ated using the Bondy-Hole-Lyttlet o n formu 



laej Bondi fc Hole I Il944t iBondi 
19531) : 



1952; McCrea 



M e = 5- 



(2GM) 2 



\3/2' 



(67) 



where S is a dimcnsionless factor of the order of 
unity. When one of the velocities, Ooo or v^, far 
exceeds the other, the accretion rate is determined 
by the dominating velocity, and can be written in 
a convenient form as (|55|) . 

In the real astrophysical situation, the param- 
eters of the surrounding matter at distances R 3> 
Rg can be taken as conditions at infinity 7 . 

As already noted, the matter surrounding a 
NS or a WD is almost always in the form of 
a high-temperature plasma with a high conduc- 
tivity. Such accreting plasma must interact effi- 
ciently with the magnetic fiel d of the compact star 
( AmnueF fc Guseinov 19681 ). Hence, the interac- 
tion between the compact star and its surround- 
ings cannot be treated as purely gravitational and 
therefore the accretion is not a purely gas dynamic 
process. In general, such interaction should be 
described by the magneto hydrodynamical equa- 
tions. This makes the already complicated pic- 
ture of interaction of the compact star with the 
surrounding medium even more complex. 

The following classification of magnetic ro- 
tators is based on the essential characteristics 



of the interaction of the plasma surrounding 
them with their electroma gnetic field. This 



(1970) 



approach was proposed by IShvartsman 
who distinguished three stages of interaction 
of magnetic rotators: the ejection stage, the 
pro peller stage, which was later rediscovered 



by llllarionov fc Sunvaev 



such, 



(|1975l ) and named as 
a nd the accret i on stag e. Using this ap- 



proach, IShvartsman 



(|l971ah was able to pre- 
dict the phenomenon of accreting X-ray pulsars 
in binary systems. New interaction regimes dis- 
covered later have led t o a general cla s sifica- 
tion of magnetic rotators ( j Lipunov Ill982al Il984 ; 



Kornilov fc Lipunov 1983aj T 



It should be noted that the interaction of the 
magnetic rotator with the surrounding plasma is 
not yet understood in detail. However, even the 
first approximation reveals a multitude of interac- 
tion models. To simplify the analysis, we assume 
the electromagnetic part of the interaction to be 
independent of the accreting flux parameters, and 
vice versa. 

Henceforth, we shall assume in almost all cases 
that the int rinsic magnetic field of a rotator is a 
dipole field (jLandau fc Lifshiz Ill97ll ): 



(68) 



This is not just a convenient mathematical sim- 
plification. We will show that the magnetoplasma 
interaction takes place at large distances from the 
surface of the magnetic rotator, where the dipole 
moment makes the main contribution. Moreover, 
the collapse of a star into a NS is known to 
"cleanse" the magnetic field. Indeed, the conser- 
vation of magnetic flux leads to a decrease of the 
ratio of the quadrupole magnetic moment q to the 
dipole moment /i in direct proportion to the radius 
of the collapsing star, q//j, oc R. 

It should be emphasized, however, that the con- 
tribution of the quadrupole component to the field 
strength at the surface remains unchanged. 

The light cylinder radius is the first important 
charactcristric of the rotating magnetic field: 



Ri = 



(69) 



7 R G 



where c is the speed of light. 

A specific property of the field of the rotating 
magnetic dipole in vacuum is the stationarity of 
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the field inside the light cylinder and formation of 
magneto dipole radiation beyond the light cylin- 
der. The lumino sity of the magnetic dipo le radia- 
tion is equal to ( Landau fc Lifshiz~lll97ll ): 



sin (3 = fci-^g lu 



(70) 



rn ~ 3 C 3 

where kt = | sin 2 /3. 

This emission exerts a corresponding braking 
torque 



K m = -|- 



2\ n 2 sin 2 (3 



(71) 



leading to a spin down of the rotator. Although 
magnetic dipole radiation from pulsars do not ex- 
ists, almost all models predict energy loss quantity 
near this value. At the same time we do not take 
into account possibly complicated angular depen- 
dence of such loss. 

9.3. The Stopping Radius 

Now we consider qualitatively the effect of the 
electromagnetic field of a magnetic rotator on the 
accreting plasma. Consider a magnetic rotator 
with a dipole magnetic moment /i, , rotational fre- 
quency lo, and mass M . At distances R Rq 
the surrounding plasma is characterized by the 
following parameters: density poo, sound veloc- 
ity doc and/or velocity Voo relative to the star. 
The plasma will tend to accrete on to the star un- 
der the action of gravitation. The electromagnetic 
field, however, will obstruct this process, and the 
accreting matter will come to a stop at a certain 
distance. 

Basically, two different cases can be considered: 

• When the interaction takes place beyond the 
light cylinder R sfn n > Ri - This case first 
considered bv lShvartsmanl (|l970l[l971al) . In 
this case the magnetic rotator generates a 
relativistic wind consisting of a flux of differ- 
ent kinds of electromagnetic waves and rel- 
ativistic particles. The form in which the 
major part of the rotational energy of the 
star is ejected is not important at this stage. 
What is important is that both relativistic 
particles and magnetic dipole radiation will 
transfer their momentum and hence exert 



pressure on the accreting plasma. Indeed, 
random magnetic fields are always present 
in the accreting plasma. The Larmor radius 
of a particle with energy <C 10 10 eV mov- 
ing in the lowest interstellar magnetic field 
~ 10~ 6 G is much smaller than the charac- 
teristic values of radius of interaction, so the 
relativistic wind will be trapped by the mag- 
netic field of the accreting plasma and thus 
will transfer its momentum to it. 

Thus, a relativistic wind can effectively im- 
pede the accretion of matter. A cavern is 
formed around the magnetic rotator, and the 
pressure of the ejected wind P m at its bound- 
ary balances the ram pressure of the accret- 
ing plasma P a : 



Pm{Rstop) — Pa{R. 



stop ) 



(72) 



This equality defines a characteristic size 
of the stopping radius, which we call the 
Shvartzman radius Rsh- 

• The accreting plasma penetrates the light 
cylinder R st0 p < Ri ■ The pressure of the ac- 
creting plasma is high enough to permit the 
plasma to enter the light cylinder. Since the 
magnetic field inside the light cylinder de- 
creases as a dipole field, the magnetic pres- 
sure is given by 



R2 ,,2 

P - _ ~ t 
m Sir ~ 8vri?6 : 



(73) 



Matching this pressure to the ram pressure 
of the accreting plasma yields the Alfven ra- 
dius Ra- 

The magnetic pressure and the pressure of 
the relativistic wind can be written in the 
following convenient form: 



Pm = 



R > R U 



We introduce a dimcnsionless factor k t such 
that the power of the ejected wind is 



(75) 
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Assuming fc t = 1/2 we get for R = Ri a 
continuous function P m — P m (R). 

The accreting pressure of plasma outside the 
capture radius is nearly constant, and hence 
gravitation does not affect the medium pa- 
rameters significantly. In contrast, at dis- 
tances inside the gravitational capture ra- 
dius Rg the matter falls almost freely and 
exerts pressure on the "wall" equal to the dy- 
namical pressure. For spherically symmetric 
accretion we obtain 



Pa = 



AnR 2 , ■ 



R > R G , 
'«fi) 1/2 , 



R<R G , 



(76) 



Here we used the continuity equation M c = 
AttRqPooVoo. When presented in this form, 
the pressure P a is a continuous function of 
distance. 

Summarizing, for the stopping radius we get 



R 



stop 



Ra, Rstop < Rl, 



Rsh, Rstop > Ri, 
The expressions for the Alfven radius are 



(77) 



Ra 



2p 2 G 2 M 2 ^ 1/6 

M C V* 



2M C V2GM 



2/7 



Ra > Rg, 
Ra < Rg, 



(78) 



and for the Shvartzman radius: 



Rsh 



, Rsh > Rg , 
(79) 



9.4. The Stopping Radius in the Super- 
critical Case 

The estimates presented above for the stopping 
radius were obtained under the assumption that 
the energy released during accretion does not ex- 
ceed the Eddington limit, so we neglected the re- 
verse action of radiation on the accretion flux pa- 
rameters. 



Now, we turn to the situation where one can- 
not neglect the radia t ion pre ssure. Consider this 



effect after iLipunov 



( 1982bl ). Suppose that the 
accretion rate of matter captured by the magnetic 
rotator is such that the luminosity at the stopping 
radius exceeds the Eddington limit (see equation 



We shall assume, following lShakura fc Sunvaev I 
([19731 ). that the radiation sweeps away exactly 
that amount of matter which is needed for the 
accretion luminosity of the remaining flux to be 
of the order of the Eddington luminosity at any 
radius: 



M(R) 



GM 



J eddi 



This yields 



(80) 



(81) 



where R s is a spherization radius (where the accre- 
tion luminosity first reaches the Eddington limit), 
and k designates the specific opacity of matter. 
Using the continuity equation, the ram pressure of 
the accreting plasma is now obtained as another 
function of the radial distance 



P a ~ P y^^v= Mc /f II R-^, R<R. 



AttR 2 



(82) 

in contrast to the subcritical regime when P a cx 
R- 5 / 2 . Matching P a and P m (see the previous 
section) for the supercritical case gives 



Ra 



= ( ^ K 



2/0 



I? 



.S7i 



\8irc^2GM x 
~ \4irc 5 V2GM J 



(83) 



(84) 



The critical accretion rate M cr is defined by the 
boundary of the inequality 



M r > M r , 



and, correspondingly, is 



4ttc 

M cr = R st . 

K 



(85) 



(86) 
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The dependence of the Alfven radius on the ac- 
cretion rate is such that the Alfven radius (be- 
yond the capture radius) slightly decreases with 
increasing accretion rate as M -1 / 6 , while it de- 
cerases below the capture radius as M~ 2 / 7 and 
attains its lowest value for the critical accretion 
rate M c > M cr , beyond which it is independent of 
the external conditions. 

We also note that in the supercritical regime, 
the pressure of the accreting plasma increases 
more slowly (as i? -3 / 2 ) when approaching the 
magnetic rotator than the pressure of the relativis- 
tic wind (as R~ 2 ) ejected by it. This means that 
in the supercritical case a cavern may exist even 
below the capture radius. 

The estimates presented here, of course, are 
most suitable for the case of disk accretion. In 
fact, the supercritical regime seems to emerge 
most frequently under these conditions. This can 
be simply understood. Indeed, the accretion rate 
is proportional to the square of the capture radius 
M c oc R G . At the same time, the angular momen- 
tum of the captured matter is also proportional to 
Rq. Hence, at high accretion rates the formation 
of the disk looks natural. 

9.5. The Effect of the Magnetic Field 

Apparently, the magnetic field of a star be- 
comes significant only when the stopping radius 
exceeds the radius of the star, R st > R x - We take 
the Alfven radius Ra for R st , since it is the small- 
est of the two quantities Ra and Rsh- Hence, we 
can estimate the lowest value of magnetic field of 
a star which will influence the flow of matter 



Mr, 



[( 



4G 2 A-/2 

(McVgm, 



, Ra > Rg, 



R 



7/2 



1/2 

j Ra < Rg 

M c > M cr , 



Mr < Mr 



(87) 



The case Ra > Rg, M c > M cr is considered 
most frequently, and for this case we get the fol- 
lowing numerical estimations 



Umin — 1U it 6 M 17 m x ' (j cm , 



or, equivalently, 

B min = 10 r Rl /A Ml( 2 m 1 J 4 G, (89) 



Most presently observed NS have magnetic 
fields ~ 10 12 G and dipole moments ~ 10 30 G 
cm , so the magnetic field must necessarily be 
taken into account when considering interaction 
of matter with these stars. 

9.6. The Corotation Radius 

The corotation radius is another important 
characteristics of a magnetic rotator. Suppose 
that an accreting plasma penetrates the light 
cylinder and is stopped by the magnetic field at a 
certain distance R s t given by the balance between 
the static magnetic field pressure and the plasma 
pressure. Suppose that the plasma is "frozen" in 
the rotator's magnetic field. This field will drag 
the plasma and force it to rotate rigidly with the 
angular velocity of the star. The matter will fall 
on to the stellar surface only if its rotational ve- 
locity is smaller than the Keplerian velocity at the 
given distance R s t'- 



wR st < y/GM x /R 3 



(90) 



1972; Davidson fc Ostriker I 



Illarionov fc Sunvaev 1 119751 ) 



Otherwise, a centrifugal barrier emerges and 
the rapidly rotati ng magnetic field impedes the ac- 
cretion of matter dShvartsman 1970; Pringle & Rccs 

^ 197.4 [Lamb I 1197.4 ~ 
The latter authors 
assumed that if LoR st \J GM X / R s t , the mag- 
netic field throws the plasma back beyond the 
capture radius. They called this effect the "pro- 
peller" regime. In fact, matter may not be shed 
([Lipunov Ill982af) . but it is important to note that 
a stationary accretion is also not possible. 
The corotation radius is thus defined as 

R c = (GM x /lo 2 ) 1/3 ~ 2.8xl0 8 my 3 (P/ls) 2/3 cm, 

(91) 

where P is the rotational period of the star. 

If R s t < R c , rotation influences the accretion 
insignificantly. Otherwise, a stationary accretion 
is not possible for R st > R c - 

9.7. Nomenclature 

The interaction of a magnetic rotator with the 
surrounding plasma to a large extent depends on 
the relation between the four characteristic radii: 
the stopping radius, R st , the gravitational capture 
radius, Rq, the light cylinder radius, Ri, and the 
corotation radius, R r . The difference between the 
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interaction regimes is so significant that the mag- 
netic rotators behave entirely differently in differ- 
ent regimes. Hence, the classification of the inter- 
action regimes may well mean the classification of 
magnetic rotators. The classification notation and 
terminology is described bel ow and summari zed in 
Table 1, based on paper bv iLipunov I (jl987t ). 

Naturally not all possible combinations of the 
characteristic radii can be realized. For example, 
the inequality Ri > R c is not possible in princi- 
ple. Furthermore, some combinations require un- 
realistically large or small parameters of magnetic 
rotators. Under the same intrinsic and external 
conditions, the same rotator may gradually pass 
through several interaction regimes. Such a pro- 
cess will be referred to as the evolution of a mag- 
netic rotator. 

We describe the classification by considering an 
idealized scenario of evolution of magnetic rota- 
tors. Suppose the parameters p^, u,^ and M c of 
the surrounding medium remain unchanged. We 
shall also assume for a while a constancy of the ro- 
tator's magnetic moment fj,. Let the potential ac- 
cretion rate M c at the beginning be not too high, 
so that the reverse effect of radiation pressure can 
be neglected, M c < M cr . We also assume that 
the star initially rotates at a high enough speed to 
provide a powerful relativistic wind. 

Ejectors (E). We shall call a magnetic rotator 
an ejecting star (or simply an ejector E) if the pres- 
sure of the electromagnetic radiation and ejected 
relativistic particles is so high that the surround- 
ing matter is swept away beyond the capture ra- 
dius or radius of the light cylinder (if Ri > Rg)- 



Ejector: Rsn > max(Ri, Rg), 



(92) 



It follows from here that P m oc R~ 2 while the 
accretion pressure within the capture radius is 
P a oc i? 5 / 2 i.e. increases more rapidly as we ap- 
proach an accreting star. Consequently, the radius 
of a stable cavern m ust exceed the capture radius 
(|Shvartsmanlll970h . 



It is worth noting that the reverse transition 
from the propeller (P) stage to the ejector (E) 
stage is non-symmetrical and occurs at a lower 
period (see below). This means that to switch a 
pulsar on is more difficult than to turn it off. This 
is due to the fact that in the case of turning-on of 
the pulsar the pressures of plasma and relativis- 



tic wind must be matched at the surface of the 
light cylinder, not at the gravitational capture ra- 
dius. In fact, the reverse transition occurs under 
the condition of equality of the Alfven radius to 
the radius of the light cylinder (Ra — Ri)- 

It should be emp hasized that, as mentioned by 
(|Shvartsmanlll970l ). relativistic particles can be 
formed also at the propeller st age by a ra pidly 
rotating magnetic field (see also Kundt I (|l990h ). 

Propellers (P). After the ejector stage, the 
propeller stage sets in under quite general condi- 
tions, when accreting matter at the Alfven sur- 
face is hampered by a rapidly rotating magnetic 
field of the magnetic rotator. In this regime the 
Alfven radius is greater than the corotation ra- 
dius, Ra > R c - A finite magnetic viscosity causes 
the angular momentum to be transferred to the 
accreting matter so that the rotator spins down. 
Until now, the propeller stage is one of the poorly 
investigated phenomena. However, it is clear that 
sooner or later the magnetic rotator is spin down 
enough for the rotational effects to be of no im- 
portance any longer, and the accretion stage sets 
in. 

Accretors (A). In the accretion stage, the 
stopping radius (Alfven radius) must be smaller 
than the corotation radius Ra < R c - This is the 
most thoroughly investigated regime of interac- 
tion of magnetic rotators with accreting plasma. 
Examples of such systems span a wide range of 
bright observational phenomena from X-ray pul- 
sars, X-ray bursters, low-mass X-ray binaries to 
most of the cataclysmic variables and X-ray tran- 
sient sources. 

Georotators (G). Imagine that the star be- 
gins rotating so slowly that it cannot impede 
the accretion of plasma, i.e. all the conditions 
mentioned in the previous paragraph are sat- 
isfied. However, matter still can not fall on 
to the rotator's surface if the Alfven radius 
i s larger than the gravitational captur e radiu s 



(jlllarionov fc Sunvaev 1975 : Lipunov 1982c ). 



This means that the attractive gravitational force 
of the star at the Alfven surface is not signifi- 
cant. A similar situation occurs in the interaction 
of solar wind with Earth's magnetosphere. The 
plasma mainly flows around the Earth's magne- 
tosphere and recedes to infinity. This analogy ex- 
plains the term "georotator" used for this stage. 
Clearly, a georotator must either have a strong 
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Table 1: Classification of neutron stars and white dwarfs. 



Abbrevi- 


Type 


Characteristic 


Accretion 


Well known 


ation 




radii relation 


rate 


observational appearances 


E 


Ejector 


R s t > Rg 
Rst > Ri 


M < M cr 


Radio pulsars 


P 


Propeller 


R c > Rst 
Rst < Rg 
R s t < Ri 


M c < M cr 


? 


A 


Accretor 


Rst < Rg 
Rst < Ri 


M c < M cr 


X-ray pulsars 


G 


Georotator 


Rg < Rst 
Rst < Rc 


M c < M cr 


? 


M 


Magnetor 


Rst > a 
R c > al 


M c < M cr 


AM Her, polars, 

soft gamma repeaters, 

anomalous X-ray pulsars 


SE 


Super- 
ejector 


Rst > Ri 


M c > M cr 


7 


SP 


Super- 
propeller 


R c < R st 
R s t < Ri 


M c > M cr 


? 


SA 


Super- 
accretor 


Rst < Rc 
Rst < Rg 


M c > M cr 


? 



magnetic field or be embedded in a strongly rar- 
efied medium. 

Magnetors (M). When a rotator enters a bi- 
nary system, it may happen that its magneto- 
sphere engulfs the s econdary star. Suc h a re gime 



was considered by iMitrofanov et al~ (1977) for 



WD in close binary systems called polars due to 
their strongly polarized emission. In the case of 
NS, magnetors M may be realized only under the 
extreme condition of very close binaries with no 
matter within the binary separation. 

Supercritical interaction regimes. So far, 
we have assumed that the luminosity at the stop- 
ping surface is lower than the Eddington limit. 
This is fully justified for G and M regimes since 
gravitation is not important for them. For types 
E, P, and especially A, however, this is not always 
true. The critical accretion rate for which the Ed- 
dington limit is achieved is 



M rr = 1.5 x 10~ 6 i? s M« 



yr 



(93) 



where R$ = R s t /10 8 cm is the stopping radius 
(Schwartzman radius or Alfven radius, see above). 
We stress here that the widely used condition 



of supercritical accretion rate 

M > 1Q- 8 (M/M q )Mq yr" 1 

is valid only for the case of non-magnetic NS, 
where R st ~ 10 km coincides with the stellar ra- 
dius. In reality, for a NS with a typical magnetic 
field of 10 11 — 10 12 G, the Alfven radius reaches 
10 7 — 10 8 cm, so much higher accretion rates are 
required for the supercritical accretion to set in. 
The electromagnetic luminosity released at the NS 
surface, however, will be restricted by L e dd, and 
most of the li berated energy may be ca rried away 
by neutrinos ( Basko fc Sunvaevl Il975f ) (see also 
section about hyper accretion). 

Most of the matter in the dynamic model of 
supercritical accretion forms an outflowing flux 
covering the magnetic rota tor by an opaque shell 
(jShakura fc Sunvaev Ill973l ). The following three 
additional types are distinguished, depending on 
the relationship between the characteristic radii: 
superejector (SE), superpropeller (SP) and super- 
accretor (SA). 
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9.8. A Universal Diagram for Gravimag- 
netic Rotators 

The classification given above was based on re- 
lations between the characteristic radii, i.e. quan- 
tities which cannot be observed directly. This 
drawback can be removed if we note that the light 
cylinder radius i?;, Shvartzman radius R$h and 
corotation radius R c are functions of the well- 
observed quantity, rotational period of the mag- 
netor p. Hence, the above classification can be 
reformulated in the form of inequalities for the ro- 
tational period of a magnetic rotator. 

One can introduce two critical periods pe and 
Pa such that their relationship with period p of a 
magnetic rotator specifies the rotator's type: 



p <p E , -> 
Pe <P < PA, 
p> PA, -> 



E or SE, 

-► P or SP, 

A, SA, G or M, 



(94) 



The values of pe and pa can be determined 
from Table 2 which defines the basic nomencla- 
ture, and are functions of the parameters Woo, M c , 
fi and M x . The parameters p and \i character- 
ize the electromagnetic interaction, while M c de- 
scribes the gravitational interaction. Instead of 
M c we introduce the potential accretion luminos- 
ity L 



L 



Rx 



(95) 



The physical sense of the potential luminosity 
is quite clear: the accreting star would be ob- 
served to have this luminosity if the matter for- 
mally falling on the gravitational capture cross- 
section were to reach its surface. 

Approximate expressions for critical periods 
(jLipunov 1119921 ) are: 



0.42^ 1/4 4^38 1/4 s, (a) 
PE={ 1.8^ 6 mV 3 /^%~ 8 1/6 b,C8) 

1.4 • lO- 2 ™- 1 /^/ 9 s , (7) 

(96) 

We assume the next indication in this formula: 

(a), M c < M cr , p < pgl ! (/?), M c < M cr , p > 
p GL - (7), M c > M cr . 



-1/4 



.38 



PA = 



1.2m- 5 /^% 7 L S8 



s, (a) 
s, 08) 
s, (7) 



3/7 



(97) 



We assume the next indication in this formula: 
(a), M c < M cr and R A > Ra] W), M c < 
M cr and R a <Rg- (7), M c > M cr . 

Here a new critical period pgl was introduced 
from the condition Rq = Rf. 



Pgl 



4kGM x 



500m x vj 2 s, (98) 



Treating the rotator's magnetic dipole moment 
\x and M x as parameters, we find that an over- 
whelming majority of the mag netor's stages ca n 
be shown on a "p-L" diagram (jLipunov Hl982aT ). 
The quantity L also proves to be convenient be- 
cause it can be observed directly at the accretion 
stage. 

9.9. The Gravimagnetic Parameter 

By expecting the expression for the stopping 
radius in the subcritical regime (M c < M cr ) one 
can note that the magnetic dipole moment \i and 
the accretion rate M c always appear in the same 
combination, 



y 



Mc 



(99) 



as was noticed by iDavies fc Pringle" ( 1981 ) . The 
parameter y characterizes the ratio between the 
gravitational and magnetic "properties" of a star 
and will, therefore, be called the gravimagnetic pa- 
rameter. Two magnetic rotators having quite dif- 
ferent magnetic fields, subjected to different ex- 
ternal conditions but with identical gravimagnetic 
parameters, have similar magnetospheres, as long 
as the accretion rate is quite low (M c < M cr ). 
Otherwise, the flux of matter near the stopping 
radius no longer depends on the accretion rate at 
a large distance. 

In fact, the number of indepe ndent parameter s 
can be further reduced (see e.g. Lipunov ( 19921) ) 
by introducing the parameter 



Y 



M c Vq 



(100) 
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Table 2: Parameter of the evolution equation of a magnetic rotators. 



Parameter 


Regime 












E, SE 


P, SP 


A 


SA 


G M 


M 








M c 


M C (R A /R S ) 


M c 


«t 


- 2/3 


< 1/3 


-1/3 


-1/3 


-1/3 -1/3 


Rt 


Ri 




R c 




i?A a 



Plotting the rotator's period p versus Y- 
parameter we can draw a somewhat less obvious 
but more general classification diagram than the 
"p-L" diagram discussed above. This permits us 
to show on a single plot the rotators with key 
parameters M c , /i and spanning a very wide 
range. 

In the case of supercritical accretion, another 
characteristic combination is found in all the ex- 
pressions: 



Y, 



(101) 



In analog to the subcritical "p-Y" diagram, a 
supercritical "p-Y s " diagram can be drawn. 

10. Evolution of Magnetic Rotators 

The evolution of a magnetic rotator, which de- 
termines its observational manifestations, involves 
the slow changing of the regimes of its interaction 
with the surrounding medium. Such an approach 
to the evolu t ion was developed in the 1970s by 
Shvartsman | (|1970) : Bisnovatvi-Kog an fc Romberg 



Illarionov & Sunyacv (11975); IShakura 



19751); IWickramasinghe fc Whelari I (|l975f ). 
Savoniie fc van den Heuvel I (|1977j ) and others. 
Three stages were mostly considered in these pa- 
pers: ejector, propeller and accretor. All these 
stages can be described by a unified evolutionary 
equation. 

10.1. The evolution equation 

Analysis of the nature of interaction of a mag- 
netized star with the surrounding plasma allows 
us to write an approximate evolution equation for 
the angular mom entum of a mag netic rotator in 
the general form ( Lipunov Ill982g 



dluj 



Ki 



RV 



(102) 



where k su is a specific angular momentum applied 
by the accretion matter to the rotator. This quan- 
tity is given by 

(GMxRd) 1 / 2 , Keplerian disk accretion, 
rjtflRQ, wind accretion in a binary, 
— 0, a single magnetic rotator. 

(103) 

where Rd is the radius of the inner disk edge, 
f2 is the rotational frequency of the binary sys - 
tem, and rjt ~ 1/4 ()lllarionov fe Sunvaev 19751) . 
The values of dimensionless factor K t , characteris- 
tic radius Rt and the accretion rate M in different 
regimes are presented in Table 2. 

The evolution equation (|102[) is approximate. 
In practice, the situation with propellers and su- 
perpropellers is not yet clear. In Table 2 R m is 
the size of a magnetosphere whose value at the 
propeller stage is not known accurately and which 
may differ significantly from the standard expres- 
sions for the Alfven radius. 

10.2. The equilibrium period 

The evolution equation presented above indi- 
cates that an accreting compact star must en- 
deavor to attain an equilibr ium state in which the 
resul t ant torque vanishes foavidson fc Ostriker 



19731 : iLipunov fc Shakura Ill976fh This hypothe- 
sis is confirmed by observations of X-ray pulsars. 

By equating the right-hand side of equation 
() 1 03(1 to zero, we obtain the equilibrium period: 



7.8^vW^(GM E )- 5 /V 3 / 7 s, (a) 

(P) 

(104) 



Peg 

p eq = V* 'A/B W L^ 7 1/2 T W 1/6 
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where A w 5 x 10 4 (3K t )^§ / 45 1 m ;r 1 s yr 1 , and 
B w S3 5.2xl0- 6 i?|mo /3 /(10 2/3 m2)/ 45 1 M_ 6 ?7 s yr" 1 , 
£ 3 7 = £/10 37 , T10 = T/10 days; (a), disk accre- 
tion; (/3), quasi-spherical accretion. 

Alternatively: 



p eq m 1.0L 37 3/7 /i3o 7 S i disk , 
p eg = 10 % 7 1/2 AC 6 1/2 x 
x(m 2/3 /(10 2 / 3 m 2 ))-i/2 x 

-1^-1/6. .__ „ stellar wind, 



(105) 



„ r — J-T-i— 1/6 
X - L 37 J 10 M30 



S, 



Let us turn to the case of disk accretion. The 
above model of the spin-up and spin-down torques 
possesses an unexpected property. The equilib- 
rium period obtained by setting the torque to zero 
is connected with the critical period pa through a 
dimensionless factor: 



.1/2 



(106) 



The parameters Kt and e must be such that 
Peq > PA- Since Kt w e ps 1, the equilibrium period 
in the case of disk accretion is close to the critical 
period, pa, separating accretion stage A and the 
propeller stage P. In the case of the supercritical 
accretion the equilibrium period is determined by 
formula 

(|Lipunovlll982bf ): 



Peq(SA) ~ 0.17a* 



2/3 _i/9 



.'ill 



(107) 



10.3. Evolutionary Tracks 



The evolution of NS in binaries must be stud- 
ied in conjunction with the evolution of nor- 
mal stars. This problem was discussed q ualita - 



tively by iBisnovatvi-Kogan fc Romberg ( 19761 ): 



Savo niie fc van den Heuvel 
Lipunov I (|l982a[ ) 



11977D, 

and other. We begin with the 
qualitative analysis presented in the latter of these 
paper. 

The most convenient method of analysis of NS 
evolution is using the "p-L" diagram. It should 
be recalled that L is just the potential accretion 
luminosity of the NS. This quantity is equal to the 
real luminosity only at the accretion stage. 

In Figure 1 we show the evolutionary tracks of 
a NS. As a rule, a NS in a binary is born when 
the companion star belongs to the main sequence 



(loop- like track). During the first 10 5 — 10 7 years, 
the NS is at the ejector stage, and usually it is 
not seen as a radiopulsar since its pulse radia- 
tion is absorbed in the stellar wind of the normal 
star. The period of the NS increases in accordance 
with the magnetic dipole losses. After this, the 
matter penetrates into the light cylinder and the 
NS passes first into the propeller stage and then 
into the accretor stage. By this time, the normal 
star leaves the main sequence and the stellar wind 
strongly increases. This results in the emergence 
of a bright X-ray pulsar. The period of the NS 
stabilizes around its equilibrium value. Finally, 
the normal star fills the Roche lobe and the accre- 
tion rate suddenly increases; the NS moves first 
to the right and then vertically downward in the 
"p-L" diagram. In other words, the NS enters the 
supercritical stage SA (superaccretor) and its spin 
period tends to a new equilibrium value (sec equa- 
tion fT07|)V 

After the mass exchange, only the helium core 
of the normal star is left (a WR star in the case 
of massive stars), the system becomes detached 
and the NS returns back to the propeller or ejec- 
tor state. Accretion is still hampered by rapid NS 
rotation. This is probably the reason underlying 
the absence of X-ray pulsars in pairs with Wolf- 



Rayet stars ( LipunovTll982cf l. Since the helium 
star evolves on a rather short time-scale (« 10 s 
yr), the NS does not have time to spindown con- 
siderably: after explosion of the normal star, the 
system can be disrupted leaving the old NS as an 
ejector, i.e. as a high-velocity radio pulsar. 

The "loop-shaped" track discussed above can 
be written in the form: 

• I+E -> I+P -> II+P -> 11+ A -> III+SA -> 
IV+P — > E+E (recycled pulsar) — > . . . 

• I+E -► I+P -> II+A -> III+SA -► IV+E 
— > (recycled ejector) — > IV+P — > E+E (re- 
cycled pulsar) — > . . . 

• Another version of the evolutionary track of 
a NS formed in the process of mass exchange 
within a binary system is: 

III+SE -> III+SP -> IV+P -> E+E ->■... 

The overall lifetime of a NS in a binary system 
depends on the lifetime of the normal star and on 
the parameters of the binary system. However, 
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Fig. 1. — Tracks of NS on the period (p) - gravimagnetic parameter (Y) diagram: track of a single NS 
(vertical line) and of a NS in a binary system (looped line). For the second track, possible observational 
appearances of the NS are indicated. 
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the number of transitions from one stage to an- 
other during the time the NS is in the binary is 
proportional to the magnetic field strength of the 
NS. 

Figure 2 demonstrates the effect of NS mag- 
netic field decay (track (a) with and (b) without 
magnetic field decay). The first track illustrates 
the common path which results in the production 
of a typical millisecond pulsar. 

10.4. Evolution of Magnetic Rotators in 
Non-circular Orbits 

So far we have considered evolution of a mag- 
netic rotator related to single rotators or those 
entering binary system with circular orbits. This 
approximation was appropriate for the gross anal- 
ysi s of binary evolutionary sc enario performed 
by iKornilov fc Lipunov I ( 1983al lbh. This approx- 
imation is further justified by the fact that the 
tidal interaction in close binaries l eads to or - 
bital circularization in a short time ( Hut I ll~98lh . 



tailed calculations, and suggests the existence of 
two different types of bina ry systems separated by 
a crit ical eccentricity, e cr (jGnusareva fc Lipunov 



1985) 



However, the more general case of a binary 
with eccentric orbit must be considered for fur- 
ther analysis. It is especially important because 
many of the currently observed X-ray pulsars, as 
well as radio p ulsars with massive c ompanions 
PSR B125 9-63 ll Johnston et all Il992h and PSR 
B0042-73 (|Kaspi et al.l Il994h . are in highly ec- 
centric orbits around massive companions. Pre- 
viously, such studies have been performed b y 
Gnusareva fc Lipunov] (|l985l) JProkhorovl (|l987l ). 

Orbital eccentricity necessarily emerges after 
the first supernova explosion and mass expulsion 
from the binary system. In massive binaries with 
long orbital periods > 10 days, the eccentricity 
may be well co nserved unt il the second episode of 
mass exchange (|Hut l ll98lh . Here, we concentrate 



on the evolutionary consequences of eccentricity. 

10.5. Mixing types of E-P-A binary sys- 
tems with non-zero orbital eccentric- 
ity 

The impact of eccentricity on the observed 
properties of X-ray pulsars has been consid- 
ered in many papers (lAmnueF fc Guseinov 1971 



Shakura fc Sunvaev ll973tlPacini fc Shapiro Il975t 
Lipunov fc Shakura 1 11976). The most important 
consequence of orbital eccentricity for the evolu- 
tion of rotators can be understood without de- 



Consider an ideal situation when a rotator en- 
ters a binary system with some eccentricity. The 
normal star (no matter how) supplies matter to 
the compact magnetized rotator. We assume that 
all the parameters of the binary system (binary 
separation, eccentricity, masses, accretion rate, 
etc.) are stationary and unchanged. Then a crit- 
ical eccentricity e cr appears such that at e > e cr 
the rotator is not able to reach the accretion state 
in principle. Let the rotator be rapid enough ini- 
tially to be at the ejector (E) state. With other 
parameters constant, the evolution of such a star 
is determined only by its spindown. The star will 
gradually spindown to such a state that when pass- 
ing close to the periastron where the density of 
the surrounding matter is higher, the pulsar will 
"choke" with plasma and pass into the propeller 
regime. Therefore, for a small part of its life the ro- 
tator will be in a mixed EP-state, being in the pro- 
peller state at periastron and at the ejector state 
close to apastron. The subsequent spindown of 
the rotator leads most probably to the propeller 
state along the entire orbit. This is due to the 
fact that the pressure of matter penetrating the 
light cylinder Ri increases faster than that caused 
by relativistic wind and radiation, as first noted 
by IShvartsman (1971a). So it proves to be much 
harder for the rotator to pass from the P state to 
the E state than from E to P state (see the follow- 
ing section). 

The rotator will spindown ultimately to some 
period, pa, at which accretion will be possible dur- 
ing the periastron passage. Accretion, in contrast, 
will lead to a spin-up of the rotator, so that it 
reaches some average equilibrium state character- 
ized by an equilibrium period p eq defined by the 
balance of accelerating and decelerating torques 
averaged over the orbital period. If the eccentric- 
ity was zero, the rotator would be in the accretion 
state all the time. By increasing the eccentricity 
and keeping the periastron separation between the 
stars unchanged, we increase the contribution of 
the decelerating torque over the orbital period and 
thus decrease p eq . At some ultimate large enough 
eccentricity e cr the equilibrium period will be less 
than the critical period pa permitting the transi- 
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Fig. 2. — The period-gravimagnetic parameter diagram for NS in binary systems, (a) with NS magnetic 
field decay (the oblique part of the track corresponds to "movement" of the accreting NS along the so-called 
"spin-up" line), (b) a typical track of a NS without field decay in a massive binary system. 



tion from the propeller state to the accretion state 
at apastron to occur. The rotational torque ap- 
plied to the rotator, averaged over orbital period, 
vanishes, and in this sense the equilibrium state is 
achieved, but the rotator periodically passes from 
the propeller state to the accretion state. 

Thus, X-ray pulsars with unreachable full-orbit 
accretion state must exist. This means that from 
the observational point of view such binaries will 
be observed as transient X-ray sources with sta- 
tionary parameters for the normal component. 

Typically, the evolutionary track of a rotator in 
an eccentric binary is 

• E -> PE -> P -> AP, e > e cr 

• E -> PE -> P -> AP -> A, e < e cr 

this may be the principal formation channel of 
transient X-ray sources. 

10.6. Ejector-propeller hysteresis 

As mentioned earlier, the transition of the ro- 
tator from the ejector state to the propeller state 
is not symmetrical. Here we consider this effect in 
more detail. In terms of our approach, we must 
study the dependence of R st on M c . To find R st , 



we must match the ram pressure of the accreting 
plasma with that caused by the relativistic wind or 
by the magnetosphere of the rotator. This depen- 
dence R st (M c ) will be substantially different for 
rapidly (Ri < Rq) and slowly (Ri > Rq) rotat- 
ing stars (see Figure 3). One can see that in the 
case of a fast rotator, an interval of M c appears 
where three different values of R s t are possible, 
the upper value R\ corresponding to the ejector 
state and the bottom value i?3 to the propeller 
state; the intermediate value R2 is unstable. This 
means that the rotator's state is not determined 
solely by the value of M c , but also depends on 
previous behavior of this value. 

Now consider a periodic changing of M c caused, 
for example, by the rotator's motion along an ec- 
centric orbit, and large enough for the rotator to 
transit from the ejector state to the propeller state 
and vice versa. Initially, the rotator is in the ejec- 
tor state. By approaching the normal star, the 
accretion rate M c increases and reaches a criti- 
cal value Mep, where the equilibrium points R\ 
(stable point corresponding to the ejector state) 
and Ri (unstable) approach Rq (upper kink), 
where they merge (see Figure 2). After that only 
one equilibrium point remains in the system, the 
stopping radius R st jumps from w Rq down to 
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Fig. 3. — Dependence of the stopping radius R s t on the modified gravimagnetic radius R y for two possible 
relations between the light cylinder radius Ri and the gravitation capture radius Rg'- Ri > Rg (left-hand 
panel) and Ri < Rg (right-hand panel; here the ejector-propeller hysteresis becomes possible). 



i?3 < Ri , and the rotator changes to the propeller 
state. 

As M c decreases further along the orbit and 
reaches the critical value Mep once again, the 
reverse transition from propeller to ejector does 
not occur. The transition only occurs when M c 
reaches another critical value, Mpe < Mep, 
where the unstable point i?2 meets the stable pro- 
peller point i?3, and the stopping radius R s t jumps 
from ~ Ri up to R\ > Rg- It should be noted 
that for fast enough rotators, a situation is pos- 
sible when the step down from the ejector state 
occurs in such a manner that the stopping ra- 
dius R s t < R c and the rotator passes directly to 
the accretion state. The reverse transition always 
passes through the propeller stage: A — > P — > 
E. In principle, transitions from the ejector state 
to supercritical st ates SP or SA are also possible 
(|Prokhorovlll987t ). In the case of slow rotators 
(Ri > Rb) )j the "E-P" hysteresis is not possible, 
and transitions between these states are symmet- 
rical. 

10.7. E-P transitions for different orbits 

Orbital motion of the rotator around the nor- 
mal companion in an eccentric binary draws a hor- 
izontal line on the " p-Y" diagram, with the begin- 
ning at a point corresponding to M c (a p ), and the 
end at a point corresponding to M c (a a ) (here a p 
and a a are the periasrton and apastron distances, 
respectively). The length of this segment is de- 
termined by the eccentricity. Since Y oc M c , the 
rotator moves along this segment from left to right 
and back as it revolves from the apastron to the 
periastron. At each successive orbital period, this 
line slowly drifts up to larger periods. The evolu- 



tion of this system is thus determined by the or- 
der the critical lines on the diagram are crossed by 
this "line" . It is seen from the " p-Y" diagram that 
the regions with and without hysteresis are sepa- 
rated by a certain value of the parameter Y = Yj.. 
Since Y oc M c oc 1/r 2 , four different situations are 
possible depending on the relationship of the bi- 
nary orbital separation a with critical va lue 



corresponding t o Yk 



Lipunov et al. I (|1996 



3b) 



see e.g. 
)■ 



Lipunov I (jl99J); 



1. a p > a cr . In this case no hysteresis oc- 
curs and transitions E-P and the reverse 
take place in symmetrical points of the orbit. 
The rotator passes the following sequence of 
stages: E — > EP — > P — > . . . Here EP means 
a mixed state of the rotator at which it is 
in the ejector state during one part of the 
orbit, and in the propeller state during the 
rest of the orbital cycle. 

2. a a > a cr > a p . In this case, the hysteresis 
occurs at the beginning of the mixed EP- 
state (state EP^), but as the rotator slows 
down the hysteresis gradually decays and 
disappears. E-P transition for this system 
is: E ->■ EP,, — > EP — > P — > . . . 

3. a a < a cr . The hysteresis is possible in prin- 
ciple, but the shape of the transition depends 
on the eccentricity. Suppose a pulsar has 
spun down so much that the first transition 
from the ejector to the propeller occurred at 
the periastron. If the eccentricity were smal 
e < e cr , (do not confuse this e cr with the 
critical eccentricity introduced in the pre- 
vious section) the reverse transition to the 
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ejector state would not occur even at the 
apastron, and the evolutionary path would 
be E -> P -> . . . 



4. If e > e c 



the track is E 



... It should be noted that just after the first 
EP transition (as well as before the last), the 
system spends a finite time in the E and P 
states at every revolution. 

The value of e cr can be expressed through the 
orbital parameters as 



1/7 



1/7 



1/7 . 1/7 > 



(108) 



To conclude, we note that the hysteresis during 
the ejector-propeller transition may be possible for 
single radio pulsars also. For example, when the 
pulsar moves through a dense cloud of interstellar 
plasma, the pulses can be absorbed. The radio 
pulsar turns on again when it comes out from the 
cloud. The hysteresis amplitude for single pulsars 
can be high enough because of small relative ve- 
locities of the interstellar gas and the pulsar, so 
that Rq > R h 

11. Summary 

The "Scenario Machine" is the numerical code 
for theoretical investigations of statistical proper- 
ties o f binary stars, i.e. this i s population synthesis 
code (jLipunov et al. 1 ll996bh . It includes the evo- 
lution of normal stars and the evolution of their 
compact remnants. This is espe cially importan t 
for studies of the neutron stars ( Lipunovl 1992). 
We always include the most important observa- 
tional discoveries and theoretical estimations into 
our code. 
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